Hacemos una tabla de los barrios de donde vienen esos usuarios y la distancia (al centro del barrio)
users_barrio_Ikea <- user_barrio_most[user_barrio_most$user %in% users_Ikea,]
agg_barrio_Ikea <- aggregate(users_barrio_Ikea$bmax,by=list(users_barrio_Ikea$bmax),length)
Juntamos estos datos con el mapa
mapb_Ikea <- mapb
mapb_Ikea@data$nusers <- agg_barrio_Ikea$x[match(mapb_Ikea@data$GEOCODIGO,agg_barrio_Ikea$Group.1)]
Pintamos un choropleth con estos datos
mapp <- mapb_Ikea[!is.na(mapb_Ikea@data$nusers),]
xx <- mapp@data$nusers
xx <- (xx-min(xx))/(max(xx)-min(xx))
m <- leaflet(data = mapp) %>% addTiles() %>%
addPolygons(fillColor= ~ rgb(xx,0,0),
weight=1,fillOpacity = 0.8,
popup=~paste(DESBDT,"<br> nusers = ",nusers,sep=""))